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We investigate, by expenments and numerical simulations, thermal activation pro- 
cesses of Josephson tunnel junctions in the presence of microwave radiation. When 
the applied signal resonates with the Josephson plasma frequency oscillations, the 
switching current may become multi-valued temperature ranges both below and above 
the the classical to quantum crossover temperature. Switching current distributions 
are obtained both experimentally and numerically at temperatures both near and far 
above the quantum crossover temperature. Plots of the switching currents traced as 
a function of the applied signal frequency show very good agreement with a sim- 
ple anharmonic theory for Josephson resonance frequency as a function of bias 
current. Throughout, experimental results and direct numerical simulations of the 
corresponding thermally driven classical Josephson junction model show very good 
agreement. 

1 Introduction 

The Josephson tunnel junction is a physical system very well studied due to 
its simplicity and nonlinearity [1]. Statistical properties of Josephson junc- 
tions have been another subject of intense investigation through, e.g., mea- 
surements of the escape statistics from the zero-voltage state, successfully 
confirming consistency with the classic Kramers model for thermal activa- 
tion from a potential well [2, 3]. Escape measurements represent a powerful 
tool for probing the nature of the underlying potential well, and applying an 
ac field to a low-temperature system has been reported to produce anoma- 
lous switching distributions with two, or more, distinct dc bias currents for 
which switching is likely. These measurements have been interpreted as a 
signature of the ac field aiding the population of multiple quantum levels in 
a junction, thereby leading to enhancement of the switching probability for 
bias currents for which the corresponding quantum levels match the energy 



2 



N. Gr0nbech- Jensen et al. 



of the microwave photons. Work performed in this direction appeared first in 
the hterature two decades ago [4, 5]. These results have significantly attracted 
interest toward Josephson junction systems as possible basic elements in the 
field of quantum coherence and quantum computing[6, 7, 8, 9, fO], and more 
recently other investigations have further indicated that the application of 
microwaves may not be the only condition under which level quantization 
can be observed in Josephson junctions [11]. 

Within the framework of this research topic we recently reported ex- 
perimental measurements conducted on a Josephson junction, operated well 
above the so-called quantum transition temperature T^^ and direct numeri- 
cal simulations of the classical pendulum model, parameterized to mimic the 
experimental device[12]. It was found that multi-peaked switching distribu- 
tions are not unique to the quantum regime (below — hujQ/2-KkB), and, in 
fact, are manifested with the same features and under the same conditions 
in the classical regime as has been previously reported for low temperature 
measurements below T*. With the present paper we wish to contribute an 
anharmonic theory that accurately captures the bias current values of the 
observed resonant peaks in the switching distributions as a function of the 
applied frequency of the microwave field. We demonstrate agreement between 
the theory, direct numerical simulations, and experimental measurements for 
direct resonances as well as harmonic and sub-harmonic resonances at tem- 
peratures both well above and near T*. 




Fig. 1. Sketch of the physical phenomenon under investigation; a driven oscillation 
energy Eac superimposed onto thermal excitations, may cause a particle to escape 
a washboard potential. 
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Figure 1 illustrates the process under investigation: in the classical one- 

degree-of-frccdom single-particle washboard potential of the Josephson junc- 
tion [4], thermal excitations (shaded in the sketch) of energy ksT and the 
energy Eac of forced oscillations due to microwave radiation, can cause the 
particle to escape from the potential well. This process can be traced by 
sweeping the current-voltage characteristics of the Josephson junction peri- 
odically. Escape from the potential well corresponds to an abrupt transition 
from the top of the Joscphson-currcnt zero- voltage state to a non-zero voltage 
state. The statistics of the switching events, in the absence of time-varying 
perturbations, have been shown to be consistent with Kramers' model [2] 
for thermal escape from a one-dimensional potential. Since the thermal equi- 
librium Kramers model does not include the effect of non-equilibrium force 
terms, the results of the switching events generated by the presence of a 
microwave radiation on a Josephson junction can be investigated, in a ther- 
mal regime, only by a direct numerical simulation of the governing equations 
(RSCJ model) [1]. 



2 Theory 

The RSCJ model reads, 

hC (P(p h dip 



Ic sin ip = + sin Udt + N{t) . (1) 



Here, is the phase difference of the quantum mechanical wave functions 
of the superconductors defining the Josephson junction, C is the magnitude 
of junction capacitance. R is the model shunting resistance, and Ic is the 
critical current, while I^c and lacSinwdt represent, respectively, the contin- 
uous and alternating bias current flowing through the junction. The term 
N{t) represents the thermal noise-current due to the resistor R given by the 
thermodynamic dissipation-fluctuation relationship [13] 

{Nit)) = (2) 
{Nit)Nit')) = 2^S{t-t'), (3) 

with T being the temperature. The symbol, 6{t — t'), is the Dirac delta func- 
tion. Current and time are usually normalized respectively to the Josephson 
critical current Ic and to Wj^^, where = \/2eIc/'hC is the Josephson plasma 
frequency. With this normalization, the coeflicient of the first-order phase 
derivative becomes the normalized dissipation a = Jiu}Q/2eRIc- It is also con- 
venient to scale the energies to the Josephson energy Ej = Ich/2e — Jc^o/27r, 
where = /i/2e = 2.07-10~^^W^6 is the flux-quantum. Thus, the set of equa- 
tions (1-3) can be expressed in normalized form as 
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(f + a(p + sm (p = r] + rjd sin Q4T + n{T) (4) 
(n(r)) = (5) 
{n{T)n{T')) = 2ae6{T - r') , (6) 

where 6 = is the normahzed temperature. The normahzed dc and ac 
currents arc rj ~ -j^ and rjd = ■y', respectively. 

For small-ampUtude oscillations around a stable (zero-voltage) energetic 
minimum we obtain the standard relationship between resonance frequency 
and bias current, 



Op = yr^ , (7) 

where we have omitted the dissipative contribution to the resonance fre- 
quency. However, this linear resonance is not directly relevant for the dy- 
namics leading to anomalous resonant switching. Looking at Figure 1 it is 
obvious that a switching event will arise from probing the anharmonic region 
of the potential near the local energetic maximum, and we therefore must 
anticipate a depression of the resonance frequency at these large amplitudes. 
In order to quantify this notion, we will adopt the following ansatz, 

(yj = (^0 + V' , (8) 

where tpo is a constant and ip represents oscillatory motion. Inserting this 
ansatz into equation (4) (for 9 = 0) yields, 

ip + sin ipo cos ip + cos (po sin ip = V + Vd sin f2dt ~ ccip ■ (9) 

Making the single-mode assumption, ip = asm{f24t + k), k being some con- 
stant phase, we obtain the following 

sm(fio = (10) 



Jo (a 




(11) 



the functions, J„, being the zero's order Bessel function of the first kind. 
Notice that Ores for a — > 0. Since multi-peaked switching distributions 

must require some switching events to happen near the resonance and others 
to happen for larger bias currents, we can estimate that the amplitude a must 
be approximately given by 

a«7r-2sin-i-^ (12) 
Jo(a) 

which represents the phase distance from the energetic minimiim to the saddle 
point. Approximating Jo (a) by its Taylor expansion, we can arrive at the 
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simple expression between the oscillation amplitude and the applied bias 
current r], 



Inserting this approximate expression into (11) gives an explicit relationship 
between the anharmonic resonance and the bias current, relevant for the bias 
current location of the anomalous secondary peak in the switching distribu- 
tion. 



3 Experiments and Simulations 

Experiments were performed on Josephson tunnel junctions fabricated ac- 
cording to classical Nb-NbAlOx-Nb procedures [14]. The samples had very 
good current-voltage characteristics and magnetic field diffraction patterns. 
The junctions were cooled in a ^He refrigerator (Oxford Instruments He- 
liox system), providing temperatures down to 360mK. Microwave radiation, 
brought to the chip-holders by a coax cable, was coupled capacitively to the 
junctions, and the junction had a maximmn critical current of Ic — 143/iA 
and a total capacitance of 6pF from which we estimate a plasma frequency 
of llIq/2tt = 42.5GHz. From this value of the plasma frequency the classical 
to quantum crossover temperature [15] = (JiuJo/2TrkB) = 'i25mK between 
classical thermal and quantum mechanical behavior can be estimated. The 
sweep rate of the continuous current Idc was I^c = 800toj4/s, and we verified 
that the experiment was being conducted in adiabatic conditions [11]. The 
junction has a Josephson energy Ej «46.4-10~^^Jin the temperature range 
from 370mK to 1.6K, and effective resistance R — 74i7. Evaluation of the dis- 
sipation parameter was based on the hysteresis of the current- voltage charac- 
teristics of the junctions[l]. We show data for two temperatures, T w 388mK 
and T = 1.6K. 

Figure 2 shows experimentally obtained results at T = 1.6K [12]. The 
lower frames of the figure displays the switching distributions in bias current 
at different microwave frequencies. The top frame shows the relationship be- 
tween the normalized current, for which the switching distributions have their 
resonant peak, and the applied frequency (normalized to the junction plasma 
frequency). Each black marker represents one of the switching distributions. 
Also shown in figure 2 is the linear resonance of equation (7), sown as a dashed 
curve, and the anharmonic resonance of equations (11) and (13), shown as 
a solid curve. The agreement between the experimental measurements and 
the anharmonic theory of the classical model is near perfect for the available 
data points, and we emphasize that the theoretical model of equations (11) 
and (13) has no fitted parameters to adjust in the comparison. Thus, the 
consistent depression of the experimental data relative to a linear resonance 
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consideration, observed in Ref. [12], should be expected and not give rise to 
re-fitting the critical current or the plasma resonance frequency. 
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Fig. 2. Experimentally obtained switching distributions, p(??), for the microwave- 
driven junction obtained for increasing values of the drive frequency. The frequency 
data points in the uppermost plot are relative to the position of the secondary peak 
in the plots. Temperature is T = 1.6K, and bias sweep rate is / = 800rnA/s. 
Dashed curve in uppermost graph represents the linear plasma resonance of (7), 
while the solid curve represents the anharmonic resonance of (11) and (13). 



Figure 3 shows experimentally obtained results at T = 388m.K, presented 
in the same manner as the data in figure 2. The agreement between the 
experimental measurements and the anharmonic theory of the classical model 
is again near perfect for the available data points. The resonance curves shown 
in figures 2 and 3 are identical, since we have not included the resonance 
dependence on the dissipation (since dissipation is very small) and since the 
measurements indicated that plasma resonance frequency and critical current 
were unchanged in the investigated temperature range. By comparing figures 
2 and 3, we notice that there seem to be no qualitative (and hardly any 
quantitative) differences between the data obtained at the two very different 
temperatures, even though the data of figure 2 is acquired at T « 5T« and 
the data in figure 3 represent T w 1.2T*. 
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Fig. 3. Same experimental situation as described in figure 2, but data is here 
acquired at T = 388mK. 

Numerical simulations of escape in a system described by equations (4)- 
(6) corresponding to the experiments with a = 0.00845, = 115.4-10"^, 4.76- 
10~^, and continuous bias sweep rate ^ = 2.1-10"^ have also been conducted 
in order to investigate the purely classical dynamics in comparison with the 
experimental measurements. The parameters have been chosen in agreement 
with the experiments discussed above. Switching distributions (each corre- 
sponding to 1,000-10,000 events), obtained for different values of the nor- 
malized drive frequency and temperature, were obtained as a function of 
the continuous bias, and secondary resonant peaks in the distribution were 
easily obtained in the classical model by adjusting the simulated microwave 
amplitude for a given frequency. 

Figure 4 shows both experimental measurements and direct numerical 
simulations of the resonant peak location as a function of applied microwave 
frequency at T = 1.6K « Experimental data are shown as box mark- 
ers and numerically obtained data are shown as circles. As in figures 2 and 
3, dashed curves represent the linear resonance (7) while the solid curves 
are generated from (11) and (13). Experimental data for the fundamental 
resonance (labeled ^2p) in the figure are the ones from figure 2. We clearly 
observe the close agreement between theory, experiment, and simulation. We 
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Fig. 4. The functional dependencies of the driving frequency upon the location 

of the secondary peak in p{ri) obtained for subharmonic and harmonic pumping. 
Circles represent numerical results and squares experimental data. Parameters are 
as given in Figure 2. Dashed curvcsrepresent the linear plasma resonance of (7), 
while the solid curve s represents the anhaxmonic resonance of (11) and (13). 



also present data for subharmonic resonances, and here too do we find very 
close agreement between simulation and experiment. The theoretical har- 
monic and subharmonic resonance curves are the ones of equations (7), (11), 
and (13) multiplied with the indicated fraction in the figure. This simple the- 
ory seems to also predict the sub-harmonic microwave induced resonant peak 
location very well. 

We finally show, in figure 5, the data similar to the ones in figure 4 taken 
at T = 388mK « 1.2T^,. Also at this temperature do wc observe very close 
agreement between experiment, simulation, and theory, amplifying the notion 
that microwave induced switching and anomalous switching distributions can 
be understood within a classical, thermal framework. 



4 Conclusion 

In conclusion, our theory, and experiments on ac-driven, thermal escape of a 

classical particle from a one-dimcnsional potential well have shown that res- 
onant coupling (harmonic or subharmonic) between the applied microwaves 
and the plasma resonance frequency provides an enhanced opportunity for es- 
cape, and we have directly observed the signatures of such microwave-induced 
escape distributions in the form of anomalous multi-peaked escape statis- 
tics at two temperatures, T = 388mK » 1.2T* and T = 1.6K w 5T*. The 
straightforward agreement between the classical hypothesis of anomalous dis- 
tributions being directly produced by ac-induced anharmonic resonances, the 



Anomalous thermal escape in Josephson systems. 



9 



2 YTTzzzzz^TTTTZ: 




0.50 0.60 0.70 0.80 0.90 1.00 

Fig. 5. The functional dependencies of the driving frequency upon the location of 
the secondary peak as described in figure 4. The temperature is here T = 388mK. 

results of numerical simulations of the classical pendulum model of a Joseph- 
son junction, and actual Josephson junction experiments indicate a consistent 
interpretation of ac-induccd anomalous multi-peaked switching distributions 
in the classical regime of Josephson junctions. 

It is noted that previous experimental work on ac-induced escape distribu- 
tions obtained at temperatures below is consistent with the observations 
presented here. Those experiments have produced ac-induced peaks in the 
observed switching distributions, and the relevant peaks are located along- 
side the expected classical plasma resonance curve, as we have also found 
here. An important observation is that the microwave-radiation frequency 
necessary for populating an excited quantum level (hud) in a quantum os- 
cillator coincides with the classical resonance frequency of the corresponding 
classical oscillator. Thus, the switching distributions obtained from classical 
and quantum mechanical oscillators may exhibit the same microwave induced 
multi-peak signatures, which in the classical interpretation is merely due to 
resonant nonlinear effects. It is evident then that multi-peaked switching dis- 
tributions are not a unique signature of quantum behavior in the ac-driven 
Josephson junction. We finally point out that similar anomalous (resonant) 
switching has been observed both experimentally [16] and theoretically [17] 
for single- fluxon behavior in long annular Josephson junctions in an external 
magnetic field. 
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